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ABSTRACT 

We derive a first order linear differential equation describing the shape of the probability density function 
of magnetic field strengths in the quiet Sun (PDF). The modeling is very schematic. It considers convective 
motions which continuously supply and withdraw magnetic structures. In addition, a magnetic amplification 
mechanism increases the field strength up to a threshold that cannot be exceeded. These three basic ingredi- 
ents provide PDFs in good agreement with the PDFs produced by realistic numerical simulations of magneto 
convection, as well as with quiet Sun PDFs inferred from observations. In particular, the distribution is approx- 
imately lognormal, and it produces an excess of magnetic fields (i.e., a hump in the distribution) right before 
the maximum field strength. The success of this simple model may indicate that only a few basic ingredients 
shape the quiet Sun PDF. Our approach provides a concise parametric representation of the PDF, as required to 
develop automatic methods of diagnostics. 

Subject headings: convection - Sun: magnetic fields - Sun: photosphere 
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1. INTRODUCTION 



Livingston & Harvevl d 19751) and ISmithsonl d 19751) 
discovered weak polarization signals in the interior of 
supergranulation cells. These magnetic signals are known 
as Inter-Network magnetic fields (IN), Intra-Network fields 
or, simply, quiet Sun fields. It has been long conjectured 
that such signals trace a hidden component of the solar 
magnetic field having m ost of the unsigned magnetic flux 
and m a gnetic energy (e.g. , IUnndlT959tlStenfldll982tlYi et alJ 
1993[ ISanchez Almeldal [19981 l2004t ISchriiver & Title! 
20031) . This sugges tion seems to be confirmed by recent 



measurements (e.g. | Stenflol 1 1982b iFaurobert-SchollI 1 1993b 
Faurobert-Scho llet all 11995b" ISanc hez Almeida & Lite: 
2000t 



2004; 



ISocas-Navarro & Sanchez Almeida! 
Truiillo Bueno ^Fall l2004t [Manso Sa inz et all 
Sanchez Almeidal l2005lh as well as by numerica l sim- 
ulations of magneto convection (e.g. ICattaneol 119991: 
20011 IStein & Nordlundl l2002t IVogled 
2005) . The signals discovered by 



Emonet & Cattaneo 
20031: IVogler et all 
Livingston & Harvev and Smithson represent the residual 



left when a magnetic field of com plex topology is observed 
with finite angular resolution ( e.g., Emonet & Cattaneoll2001l : 
ISanchez Almeida et all [2003). Such complex field is to 
be expected as a result of the MHD interaction between 
magnetic fields and the random motions associated with the 
granulation. 

Characterizing these IN fields is therefore important, and 
significant advances have been produced during the last years 
(see the references given above). Among the observational 
parameters used for characterization, the probability den- 
sity function of magnetic field strengths (PDF) turns out 
to be particularly useful. It is defined as the fraction of 
quiet Sun occupied by magnetic fields of each strength. It 
condenses basic physical information like the filling fac- 
tor, the unsigned flux, and the magnetic energy corre spond- 
ing to each field strength (see ISanchez Almeida 2004)). In 
addition, so me of the diagnostic te c hniques employ them 
directly (e.g ., Faurobert-Scholl 19 931: iFaurobert-Scholl et al] 
11995b iFaurobert et alJl200U iTrujillo Bueno et al.ll2004l) . and 
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PDFs are predicted by the numerical simulations, allowing 
a direct comparison between simulations and observations 
dCattanedll999tlStein & Nordlundll2002tlVogler & Schussleri 
l2003UVogler et al.ll2005l) . In spite of these advantages, mea- 
suring the quiet Sun PDF is not trivial. All direct estimates 
are strongly biased. The Zeeman signals originate only in re- 
gions of significant field strength and net polarity, whereas 
the Hanle signa ls are insensitive to the fields stron ger than a 
few hundred G. Dommguez Cerde fiaet all (l2006al) carry out 
a first attempt to provide the full PDF from G to 1800 G. 
They remove the known biases by using Hanle effect mea- 
surements to constrain the shape at weak fields (< 200 G), 
Zeeman signals for the strong fields (say, > 400 G), and as- 
suming the PDF to be continuous in between. They work out 
a set of PDFs compatible with the observations and consis- 
tent with numerical simulations of magneto-convection. As 
it is acknowledged by the authors, they carry out an ex- 
ploratory estimate and therefore further independent work 
is required for confirmation. Support would be provided if 
some of the characteristic and unexpected properties of these 
PDFs can be explained in terms of simple physical mecha- 
nisms. In particular, the empirical PDFs show an increase 
of strong kG fields right before the maximum possible field 
strength, set by the gas pressure of the quiet photosphere 
(see the solid line in Fig. [Ttt with a hump at B ~ 1700 G). 
Despite the fact that these strong fields occupy only a small 
fraction of the surface, they contribute with a large part 
of the quiet Sun unsigned magnetic flux and magnetic en- 
ergy, which makes them particularly important. Moreover 
the kG fields turn out to be ideal for studying the quiet Sun 
magnetism since they show u p in unpolarized light images 
( Sanch ez Almeida et a 1.1 12004), The possible physical origin 
of such a hump is addressed here. We discuss that a ter- 
minal hump is to be expected if a magnetic amplification 
mechanism operates in the quiet Sun. These mechanisms 
have a long tradition to explain the pr esence of kG mag- 
netic fields in plage and network regions (IWeisslfr9 66: Parker 



1978 



2001 



Spruit 1979; van Ballegooiien 1984; Sanchez Almeida 



Cameron & Gal loway 2005). They concentrate weakly 



magnetized plasma conserving the magnetic flux, which de- 
mands increasing the field strength. No matter the details of 
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FIG. 1 . — (a) Set of quiet Sun PDFs derived from numerical simulations of 
magneto convection (the discontinuous lines), and a semi-empirical PDF that 
fits observed Zeeman and Hanle signals (the solid line), (b) ks repr oducing 
the PDFs shown above. The labels in the inset stand for: Vogler = [Vogler 
12003), and DC et al. = IDommguez Cerdeiia et al]<2006al) . 



the mechanism, all of them abruptly die out at the maximum 
field strength imposed by the gas pressure of the quiet photo- 
sphere. The magnetic structures tend to pile up at this limit 
giving rise to a hum p (bottleneck effect). The idea wa s al- 
ready put forward by Dommg uez Cerdena et al.l (l2006bl) . but 
it is elaborated in this paper. 

We write down a differential equation to describe the time 
evolution of the IN PDF (§ |2}. The stationary solutions of 
such equation are analyzed in § [3] It is shown how a terminal 
hump appears in a natural way. The efficiency of the amplifi- 
cation mechanism is parameterized using a coefficient (speed) 
with units of magnetic field strength per unit time. Given 
a PDF, one can infer by inversion the speed of the concen- 
tration mechanism. The technique is applied to PDFs com- 
ing from numerical simulations, as well as to semi-empirical 
PDFs (§H). Other uses and limitations of the model PDFs are 
discussed in §|5] 

2. DIFFERENTIAL EQUATION 

We use a heuristic approach to derive an equation for the 
shape of the PDF P(B). The different mechanisms creating, 
destroying and modifying magnetic fields determine P(B). 
Here we consider three of them plus the buffeting of the gran- 
ulation, namely, (1) a magnetic amplification mechanism that 
tends to increase the magnetic field strength, (2) the sub- 
mergence of the existing flux transported by the granulation 
downdrafts, (3) the emergence of new magnetic flux trans- 
ported by the granulation upflows, and (4) the buffeting of the 
granular motions on the magnetic structures. In the stationary 



state they are perfectly balanced, therefore, 
dP(B) 



dt 



dPi(B) + dP 2 (B) + dP 3 (B) + dP 4 (B) _ Q 



dt 



dt 



dt 



dt 



(1) 



where dPj/dt represents the variation with time produced by 
the y'-th mechanism. We proceed by working out the contri- 
bution of the four terms separately, and then considering all 
them together. Simplicity is the main driver therefore a num- 
ber of simplifying hypotheses restrict the problem to make it 
tractable. 

First, consider the magnetic amplification mechanism. We 
start off by assuming that the photospheric magnetic fields are 
made of fluxtubes, each one having a single magnetic field 
vector and a single cross-section. Four parameters charac- 
terize each fluxtube; the magnetic field strength B,, the mag- 
netic field inclination (?,, the magnetic field azimuth fa, and 
the section A,. Then the properties of the photospheric mag- 
netic fields are fully characterized given the four parameters 
Bi, 6j, (pi, and A, corresponding to each fluxtube - the subscript 
i varies from one to the number of fluxtubes present in the 
photosphere. By definition, and when AB — > 0, P{B)AB is 
the fraction of photospheric volume occupied by all magnetic 
structures with strengths between B— AB/2 and B+ AB/2. It 
can be expressed as a sum that considers all fluxtubes, 



P(B)AB = 



1 



$>n( 



Bi-B 
AB 



(2) 



with Vj the volume of the i— th fluxtube, and II the rectangle 
function, 

nw = {i N<iA (3 ) 

1 elsewhere. 
Given a field strength B, the sum in the numerator of equa- 
tion d2J differs from zero only when B - AB /2 < B, < B + 
AB/2 and, therefore, it gives the total volume occupied by 
fluxtubes with field strengths between B - AB / 2 and B + 
AB/2. Note also that the normalization constant in the de- 
nominator of equation (|2]) is the volume of the photosphere 
and, therefore, a constant. In the case of interest, when 
AB — * 0, the rectangle function can be replaced with a Dirac 
(5-function since, 



S(x) = lim 

Ax^O 



Tl(x/Ax)^ 
Ax ' 



(4) 



see, e.g., iBraceweM d 1978b . This replacement leads to the 
compact expression for P(B) adequate for mathematical ma- 



nipulations, 



2-i i v ' : 



(5) 



Each fluxtube has a single inclination, therefore, we are im- 
plicitly assuming that the fluxtubes are straight and longer 
than the range of heights of the photosphere; see Figure [2] 
Accordingly, the volume of photosphere occupied by a flux- 
tube is simply set by its section and inclination, 



v; = AzAi/ cos9j, 



(6) 



where Az stands for the vertical extent of the photosphere or, 
if this is too large for the approximation © to hold, Az rep- 
resents the range of heights to be described by our PDF. We 
characterize the magnetic amplification mechanism with the 
rate n at which it changes the magnetic field strength, 



dB > 

n(Bi)= —. 
dt 



(7) 
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The rate has units of magnetic field strength per unit time, and 
we will call it speed. According to equation (0, the speed 
only depends on the field strength of the fluxtube. This is 
a working hypothesis which may not hold in some practical 
cases, as we point out in §|5] All the magnetic amplification 
mechanisms proposed so far conserve the magnetic flux so 
that an increase of field strength comes together with a de- 
crease of the section of the magnetic structure to maintain the 
productA,B, constant 1 , 

d -(B i A i ) = 0. (8) 



dt 

Equations © and (O lead to, 

dvi 



(9) 



dt Bi 

where we assume that the concentration mechanism does not 
change the field inclination in a preferred sense. Using equa- 
tions (0 and (O, and the chain rule, 

dPi(B) 



dt 



E 



^SiBt-B^KBM 6 -^ 



B) 



B 



(10) 



where we have e mployed the expr ession for the derivative of 
a 5-function (e.g.. lBracewelfl978l) . 

dS(x) 



dx 



= —S(x)/x. 



Keeping in mind the following properties, 

-™ J(S ,- S) A S ,- S) , 



and 



K(Bd 



B, 



B;-B 



B 



dn(B) 



B) 



S(Bj-B) 
Bj-B ' 



(11) 



(12) 



(13) 



one can rewrite equation (fT3l as, 



S(Bi-B) 



k(B) 



dS(B r 



B) + d -^m-B), 



Bj-B " v ' dB dB 
and, consequently, equation ( fTOb turns out to yield either, 



(14) 



dPi(B) 



dt 

or, alternatively. 



, dP(B) dn(B) 

B~ l k{B)P(B) + k(B) — — + P(B) - 



dP x (B) 
dt 



dB 

1 d[K(B)BP(B)] 
B dB ' 



dB 



(15) 



(16) 



The previous equation describes the variation with time of the 
PDF due to a magnetic amplification mechanism with speed 
k{B). Note that if rather than an amplification mechanism 
one deals with a magnetic de-amplification mechanism, equa- 
tion ( [Tol l remains valid with k(B) < 0. Consequently, one can 
think of k(B) as the average between all the amplification and 
de-amplification mechanisms operating in the quiet Sun. The 
existence of de-amplifications is indeed likely since they of- 
ten come together with the amplificatio n mechanisms (e.g., 
lSpruidll979tTcameron & Gallowavll2005h . 

1 It follows from the invariance of the magnetic fl ux under stretching 
of flux tubes in perfectly conducting plasmas. See, e.g., Childress & Gilbert 
(1995, §1.2.3). 




FIG. 2. — Schematic with a portion of the photosphere (the box) which 
includes the full range of photospheric heights Az. The two cylinders repre- 
sent magnetic fluxtubes. The photospheric volume occupied by one of such 
fluxtubes depends on its cross-section A and its inclination 9, but not on its 
length if it is longer than Az/ cos 6. 

Consider the magnetic fields transported by the downdrafts. 
In this case magnetic structures disappear in proportion to the 
existing magnetic structures, i.e., 



dP 2 (B) 
dt 



-(3(B)P(B), 



(17) 



where (3(B) provides the rate of submergence per unit time. 
The rate /3(B) depends on B to consider that different photo- 
spheric field strengths might submerge at different rates. 

The emergence of magnetic fields transported by the up- 
flows has the same structure of the submergence but consider- 
ing a PDF M(B) not necessarily the same as the photospheric 
PDF. In this case, 



dP 3 (B) 
dt 



■ a(B)M(B), 



(18) 



with the rate of emergence a(B) also depending on the field 
strength. 

The interaction between granular motions and magnetic 
fluxtubes changes magnetic field inclinations, azimuths, field 
strengths, and sections. The changes of field strength and sec- 
tion are important, but they have been included in the for- 
malism among the magnetic amplification de-amplification 
mechanisms. We are left with changes of magnetic field di- 
rection, in particular, with changes of magnetic field incli- 
nation since the PDF is immune to the azimuths; see equa- 
tions © and ©. We assume, however, that the random granu- 
lar motions do not change inclinations in a preferred direction 
which, together with equations © and (O, lead to, 



dP 4 (B) 
dt 



.0. 



(19) 



Some limitations of this assumption are pointed out in § [5] 

Inserting equations ( fToT ), (1171 . ( fT8l , and (T% into equation 
(Q}, one ends up with the following differential equation for 
the shape of P(B), 



d[K(B)BP(B)] 



dB 



-f3(B)BP(B) + a(B)BM(B). (20) 



By definition, the three coefficients that characterize the three 
physical mechanisms modifying P(B) are all positive, i.e., 
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a(B) > 0, (3(B) > and k(B) > 0. As we point out above, k(B) 
stands for the speed of the amplification mechanism whereas 
a(B) and /3(B) represent inverse of time scales for the vertical 
transport of magnetic structures. 

From now on we drop from the equations the dependence 
of all variables on the field strength, which provides manage- 
able expressions without compromising clarity. Then equa- 
tion ( f20b becomes, 



d(nBP) 
dB 



-(3BP+aBM. 



(21) 



It is a first order linear differential equation and, consequently, 
its solutions are given by, 

P=(B K y l aB'Mexp[- / {(3 / n)dB"}dB' , (22) 
Jo JB' 

where we have considered that the product kBP — > when 
B — ► 0. Note that the two coefficients a and (3 cannot be in- 
dependent since P is a PDF and therefore it must be properly 
normalized, 



PdB=\. 



According to equation d22"l i. P scales with a, therefore, given 
k and (3, a scaling factor applied to a guarantees the proper 
normalization. It follows from the solution of the differential 
equation, 

d[ K Bf] 



dB 



-pBf+a f BM, 



(24) 



which is identical to the original equation d2T1 > except that a 
has been replaced with an unconstrained a./. One can readily 
show by direct substitution that, 

f °° 

P = f/ fdB, (25) 
Jo 

is a properly normalized solution of equation (TlTT i with, 

a = a f / / fdB. (26) 
Jo 

3. PROPERTIES OF THE SOLUTIONS 

3.1. Hump at large field strength 

First and most important in the context of this paper, equa- 
tion (f2Tb predicts PDFs with a hump at large magnetic field 
strengths. The only condition is an abrupt cease of the mag- 
netic amplification mechanism when B approaches a maxi- 
mum value B max . Explicitly, 



when B —> B„ 



with 



I dn I 



large enough. 



(27) 



(28) 



Equation d27l ) puts forward a very natural condition to be sat- 
isfied by any amplification mechanism. The magnetic struc- 
tures must be in mechanical balance within the photosphere 
and, therefore, the magnetic field strength cannot provide a 
magnetic pressure exceeding the (gas) pressure of the quiet 
photosphere (e.g., ISpruitl Il98 lh . At the base of the pho- 
tosphere and in th e standard ID model atmospheres (e.g. 
Maltbvetal.lU986h . 



1800 G. 



(29) 



In order to prove that conditions d27T i and (l28l i necessarily pro 
duce a hump, one rewrites equation (fJU as, 

dP ^ rf3+dn/dB 
d~B~~t k 



■ + B~ l ]P, 



(30) 



where we have assumed that the strong fields are produced by 
the amplification mechanism rather that transported upward 
by the granular flows (M ~ when B — > B max )- Consider the 
condition d27i >. It forces 

P^0 when B — ► B mRX , (31) 



since there are no sources (M 
cannot be intensified to B = B m 
and the constraint P > 0, 



; 0) and the magnetic fields 
4 . Because of the limit OTT l 



dP 

— < when B 
dB 



B n 



(32) 



Consequently, the condition for a continuous P to have a max- 
imum (dP/dB = 0) is equivalent to 



dP 

dB >0 > 



(33) 



(23) at a B < B max - Using equations fl30"l > and ( 1331 ), there is a maxi- 



mum if 



iii>[^i- 



(34) 



(35) 



These two equations provide specific expressions for what 
large enough means in equation ((28). 

3.2. PDF for small magnetic field strengths 

We expect the direction of the magnetic field vector to be 
random when B — > 0. In this case the magnetic forces can- 
not back-react on the flows and the magnetized plasma is 
freely dragged, bended and moved around by the granular mo- 
tions. This continuous buffeting impinges a random tilt to the 
magnetic fields, which should sho w no preferred orientation. 
Dommguez Cerdena et al. (2006a) argue that this random ori- 
entation forces, 



when B -> 0. 



(36) 



For a point in the atmosphere to have B = 0, the three Carte- 
sian components of the magnetic field vector must be zero 
simultaneously, and this is a very improbable event when the 
orientation of the field is random and therefore the three com- 
ponents independent. Here we go a step further and conjec- 
ture that P should follow a Maxwellian distribution, which is 
the PDF characteristic of the modulus of a vector field with 
rando m orientation (see, e.g.. lProkhorov|[T990T: lAguilar Peris! 
fl98lh . Explicitly, 



' 7T 0~ 



-B 2 exp(-B 2 /a 2 ), 



(37) 



when B — > 0. In order to obtain this particular shape, M in 
equation (l2TT > must follow a Rayleigh distribution, 



M~^Bexp(-B 2 /a 2 m ), 



(38) 



with a = <r m y/5/7> and a = l&K/5<jy/ir. This can be checked 
by direct substitution of M and P into equation (l2Tl i. consid- 
ering that k and a do not vary with the magnetic field when 
B — > 0, and neglecting the influence of the magnetic flux trans- 
ported downward (J3P <C olM when B — > 0). 
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3.3. Analytical approximation 

It was found by [Dommguez Cerden a~et alj (12006a) that 
the PDFs produced by numerical simulations of magneto- 
convection follow a lognormal distribution, and this shape 
was adopted to describe P for sub-kG magnetic field 
strengths. In this case, 



1 



ZttsB 



exp{- 



MB/B s h 



(39) 



with the parameters B s and s related to the mean and the vari- 
ance of the distribution. It turns out that equation d2"TT > predicts 
a lognormal distribution when (a) the influence of M is neg- 
ligible (e.g., when B is larger than the typical field strengths 
supplied by the granular upflows), and (b) the various coeffi- 
cients characterizing the differential equation have a weak de- 
pendence on the field strength. Equation (|2"TT i can be rewritten 
as, 

d MBP )= p djn. 

1 k dlnB i 



dlnB 



When the variation of the right-hand-side term of this equa- 
tion with magnetic field is not is not very large, then it can be 
approximated as, 



8 d\nn 
-B+- 



a + b\nB, 



(41) 



K d\nB 

with a and b two constants 2 . In this case the integration of 
equation (|40b automatically gives a lognormal (|39l ) with, 

\nB s = -a/b, 

s = y/IJb. (42) 

Therefore, provided that the approximation (l4TT i is valid, a 
lognormal yields a good analytical representation of the PDFs 
proposed in the paper. 

3.4. Magnetic flux and magnetic energy 

The two terms on the right-hand-side of equation (l2TT i rep- 
resent the rate of magnetic flux emergence and submergence 
(§ ID- Since the magnetic amplification mechanisms do not 
change the flux (§ |2}, they have to be equal once integrated 
over the full range of field strengths. This property is auto- 
matically fulfilled by the solutions of equation (f2Tb . Integrat- 
ing (I2TI 1 from B = to B = Z? max and considering the limit (|3TT >. 
then, 

8BPdB = / aBMdB, (43) 



as demanded by the conservation of magnetic flux. The 
amplification mechanisms transform weak fields into strong 
fields. The energy of the distribution scales with the square of 
the magnetic field strength, therefore, the magnetic amplifica- 
tion mechanism is expected to create magnetic energy. One 
can show that any amplification mechanism increases the en- 
ergy coming from below multiplying equation (|2TT > by B an 
integrating from B = to B = B mdx . A trivial manipulation 
leads to, 



8B 2 PdB-- 



aB MdB + 



nBPdB. (44) 



2 Consider a field strength B\ typical of the range of B where the approx- 
imation holds. Since M has to be negligible, B{ is well above zero and B can 
be expanded as a polynomial of ln(B/B[), namely, B/B\ = exp[ln(B/Bi)l = 
l + ln(B/B|) + [ln(B/Bi)] 2 /2 + ... By construction B/Bj ~ 1 and, therefore, 
the second and higher order terms of the expansion can be n eglected, which 
justifies using a logarithm to approximate B in equation )4U . 
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FIG. 3. — (a) PDFs predicted by equation )2U depending on the speed 
of the amplification mechanism (k). (b) Speeds responsible for the above 
PDFs. The legend in (a) indicates the only parameter that differs between the 
reference (the solid line) and the other PDFs. 

Keeping in mind that the second term of the right-hand-side 
is always positive (k, B and P are positive), then 



8B 2 PdB > 



aB MdB. 



(45) 



This inequality implies that magnetic energy transported 
downward (proportional to the left-hand-side term) is always 
larger than the energy coming from below (the right-hand-side 
term). 

3.5. Examples of PDFs 

In order to illustrate some of the properties described above, 
we integrate equation d2lT > under simple assumptions that still 
provide a hump at large field strengths. We consider a and 8 
to be constant, whereas k is almost constant up to a magnetic 
field B K where it abruptly drops to zero in an interval AB K , 
namely, 

..?[!— h(t^)]. (46) 

The symbol kq stands for n when B -C B^- AB^. The dif- 
ferential equation is integrated using a standard fourth-order 
Runge-Kutta method, from low to strong fields, and using as 
boundary condition that given in equation d36T ). The feeding 
fun ction M is taken to be a Rayleigh distribution as discussed 
in § 13.21 The solid line in Figure^ shows P for B K = 1600 G, 
AB K =100 G, kq/8 =500 G and a m = 50 G, and it is used as a 
reference to explore the dependence of P on k. The curves in 
Figure^ have been computed using the speeds shown in Fig- 
ure^. First, note how all PDFs show a terminal hump. The 
hump rises as the jump of k sharpens (see the dashed line, 
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with AB K =25 G), and its position moves according to B K 
(see the dotted-dashed line, with B K =1000 G). The faster the 
magnetic amplification mechanism the larger the tail of strong 
fields. The dotted line corresponds to a kq almost twice the 
reference value. Finally, the tail of strong fields also increases 
with the magnetic flux supplied by the granulation upflows. 
The triple dot-dash line is fed with three times more magnetic 
flux than the reference PDF. This may explain why plage re- 
gions have more kG fields than the quiet Sun. 

4. INVERSION OF SELECTED PDFS 

Knowing P it is possible to derive by inversion the function 
k characterizing the concentration mechanism. This is par- 
ticularly simple for the range of large field strengths, where 
one expects M to be negligible. Under this condition, equa- 
tion (f2TT > admits a solution of the kind, 



BP: 



B'PdB'- 



B'PdB'. 



(47) 



The symbol Bq stands for the smallest field strength where M 
can be neglected (M ~ for B > Bo). The coefficient (3 has 
been assumed to be constant for convenience, but this con- 
dition can be relaxed if required. Equation d47b provides k 
given P. It has been applied to a set of PDFs to illustrate 
the feasibility of the inversion. Theses PDFs are shown in 
Figure QJ, whereas the associated k are in Figure [TJ). The 
solid line corresponds to a semi-empirical PDF that repro- 
duces the l evel of Zeeman and Hanle signals observed in the 
quiet Sun ( jDom mguez Cerde na et al.l l2006aL the reference 
PDF). Despite all the uncertainties, it represents the kind of 
the PDF to be expected in the quiet Sun. The inversion pro- 
vides k//3, i.e., the speed of the concentration mechanism 
times the time-scale for the granular downdrafts to remove 
all magnetic structures from observable layers. In order to 
reproduce the semi-empirical PDF, k//3 has to be of the or- 
der of 1000 G (between 500 G and 1700 G, according to the 
solid line in Fig.QJi). In other words, during a time scale (3~ l , 
the concentration mechanism must increase the field strength 
from zero to lkG If we associate /3" 1 with the granulation 
turn-over time, the concentration mechanism must operate 
with a time-scale of, say, 15 min. This value is not very 
different from the time-scales characteristic of the concen- 
tration mechanisms proposed in the l iterature. For ex ample, 
the convective collapse simulations by Takeuchi (1999) create 
k G fields in minutes, where as the thermal relaxation proposed 
by lSanchez Alm eida (2001) easily concentrates low magnetic 
flux structures within a granulation turn-over time. The other 
PDFs in Figure QJ cor respond to nu merical simulations of 
magneto convection by Voglerj (120031) having three different 
levels of unsigned magnetic flux (see the inset). The values of 
k/P are similar to that of the semi-empirical PDF, however, 
the variation with the field strength is significantly different. 
They peak at hG field strengths whereas the semi-empirical 
k has its peak at kG field strengths. The main effect is a re- 
duction of the kG fields of the numerical PDFs with respect 



to the semi-empirical PDF. As a by product, the gradient of k 
with B is smaller i n the numerical PDF which, according to 
the discussion in § 13.11 necessarily produces humps that are 
smaller. 

5. DISCUSSION AND CONCLUSIONS 

We derive a first order linear differential equation describ- 
ing the shape of the probability density function of magnetic 
field strengths in the quiet Sun (PDF). The modeling considers 
convective motions which continuously supply and remove 
magnetic structures from the observable layers. In addition, a 
generic magnetic amplification mechanism tends to increase 
the field strength up to a maximum threshold. The basic 
ingredients that we consider yield PDFs in good agreement 
with the PDFs produced by realistic numerical simulations of 
magneto convection, as well as with quiet Sun PDFs inferred 
from observations. In particular, they give a natural expla- 
na tion for the hump right before the c utoff at 1700 G found 
by lDomfnguez Cerdena et all (2006b a). It is produced by the 
abrupt end of the amplification mechanism at the cutoff. The 
solutions of our equation provide a parametric family of PDFs 
as r equired for diagnostic s in automatic inversion procedures 
(see lSocas-Navarrol2001l) . For example, the PDFs analyzed in 
§|4]depend on four free parameters which, in principle, can be 
fitted using four independent observables. We also derive an 
approximate expression to infer the speed of the amplification 
mechanism once the PDFs are known. In order to reproduce 
the semi-empirical PDFs, the amplification mechanism must 
concentrate field strengths from zero to 1 kG in a time-scale 
similar to the granulation turn-over time. 

The magnetic structure of the quiet photosphere is certainly 
much more complex than the modeling carried out in the pa- 
per, as it is evidenced by the existing numerical simulations of 
magneto convection (see § [T). To mention just a few poten- 
tially important ingredients that we have ignored, the speed 
of the amplific ation may dep end on the section of the mag- 
netic structure (Spruitl [l979l) . we expect a tendency for the 
plasm a with large fie ld strength to be buoyant and so vertical 
(e.g. ISchiisslerl [l986). and the fluxtubes forming loops fully 
embedded in the photosphere do not change their volume un- 
der magnetic amplification. This caveat notwithstanding, one 
should not underestimate the ability of the model to repro- 
duce realistic PDFs. It is a non-trivial property which may 
reflect that, even if sketchy and primitive, the model contains 
the basic physical ingredients shaping the PDF of the quiet 
Sun magnetic fields. Obviously, we cannot discard that the 
agreement is due to a fortunate coincidence, and so the model 
is offered here only as a mere possibility to interpret the quiet 
Sun PDF. 
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